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Diesel engine modeling draws the greater attention due to its higher efficiency as compared to spark 
ignition (SI) engine. Still it is very challenging job to model diesel engine due to its complex combustion 
phenomena. The focus of the present study is to review the different available model used for modeling 
of CI engines. The modeling of CI engine is divided into single zone, multizone and multi-dimensional 
model. Which further subdivided in many submodel i.e. heat transfer, ignition delay period, droplet 
evaporation, intake and exhaust flow, chemical kinematics and soot formation model. A comparative 
study has also been carried out with experimental validation to show the compatibility with different 
modeling approach. Some optimum mathematical input parameter has been suggested by the analysis 
of different modeling approach to minimize the NO, emission and soot formation to make diesel engine 
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1. Introduction 


One of the major pollution contributors in today’s environ- 
ment is the internal combustion (IC) engine which has been 
either spark or compression ignition (CI) engine. Due to scarcity of 
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the crude oil and its high rising price makes a challenging job 
for the engine manufacturers to manufacture or model less 
polluting and more efficient engine in order to meet its future 
demands [1,2]. 

Engine simulation has been extensively used to improve the 
engine performance. Compression ignition direct injection (CIDI) 
diesel engines have been widely used in heavy-duty vehicle, 
marine transportation and now have been increasingly being 
used in light duty vehicles, particularly in Europe and Japan [3]. 
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Nomenclature 

A pre-exponential factor 

An minimum area of nozzle (m) 

Ac cylinder area (m°) 

Ae effective flow area (m?) 

ay empirical constant 

Ca coefficient of discharge 

D cylinder diameter (m) 

Ea activation energy 

Fs geometric circulation area of inlet (m?) 

K adiabatic exponent of working fluid in a cylinder 
ki rate constant (cm?/mol-s) 

Ks adiabatic exponent of working fluid before inlet valve 
m mass flow rate (kg/s) 

Ms mass of soot 

Mfg masses of vaporized fuel 

msp rates of soot formation 

P cylinder pressure (bar) 

P; preparation rate (kg/Crank angle) 

Pa pressure at the inlet (bar) 

Pm motored pressure (bar) 

Pi upstream and stagnation pressures (bar) 
P2 downstream stagnation pressures (bar) 

Pm average mean pressure 

Poz partial pressure of oxygen 

Ps pressure of working fluid before inlet valve 
P; pressure of working fluid after outlet valve 


R universal gas constant (J mol~! K7!) 
Re Reynolds number 

T temperature inside the cylinder (K) 
Tw temperature of the cylinder wall (K) 
Ta activation temperature (K) 

Tg temperature of gas (K) 

up piston speed (m/s) 

V volume of the cylinder (m?) 

Greek 

Ap pressure drop (bar) 

do change in crank angle (degree) 

oO surface tension (N/m) 

Ls inlet valve flow coefficient 

Pa density of liquid (kg/m?) 

0 spray angle (rad) 

y ratio of the specific heats 

Ag gas thermal conductivity 

la viscosity of the ambient gas 
Subscripts 

inj injection of fuel 

soc start of combustion 

f fuel 


Experimental work which is aimed at fuel economy and low 
pollutants emission for IC engine requires change in input para- 
meter which is highly demanding in terms of money and time. 
So, in order to overcome this drawback, an alternative simulation 
of engine performance with the help of mathematical model and 
powerful digital computers lowers the cost and time. In these 
simulation models, the effect of various design structures like 
design of combustion chamber input parameters (intake pressure, 
injection timing, etc.) and operation changes (compression ratio, 
speed, etc.) can be estimated in fast and non-expensive way 
provided that main mechanism are recognized and modeled 
perfectly to meet the experimental results [4,5]. 

In case of single zone model cylinder temperature, pressure 
and mass can be obtained from ordinary differential equations. 
These models do not take into account the air entrainment, 
vaporization of fuel droplet and spatial variation of mixture 
composition and temperature. Single zone models are advanta- 
geous also because of their simplicity and their wide usage along 
with empirical data within the engine industry to make design 
decisions. The measured pressure rise in an engine is used to tune 
the model and helps in calculating the rate of heat release from 
the engine cylinder [6]. 

For the accurate prediction and study the behavior of exhaust 
emission pollutants, forced the researchers to develop two 
zone combustion models [7-10]. It provides more accurate result 
as compared to single zone. Multizone model overcomes the 
drawbacks of single zone model taking into account both the 
spatial as well as temporal variation of temperature and concen- 
tration, where the detailed analysis of fuel-air distribution 
which permits calculation of exhaust gas composition within 
the reasonable accuracy [11-14]. Multizone models separate 
the fuel spray into a large, finite number of zones. Zones are 
small packets of fuel that move through the combustion 
chamber. A second type of multizone model follows the 


development of equivalence ratio zones surrounding a central 
liquid core. 

In multi-dimensional model, the instantaneous conservation 
time average equation of mass, momentum, energy has been 
taken into calculation. Stochastic as well as computational fluid 
dynamics (CFD) model are introduced to see the effect of gas 
phase turbulence on the liquid droplets. However, the effect of 
droplets on turbulence in many parts has also been ignored. 
Cylinder is considered to be divided into number of zones and 
behavior of each individual zone being studied. 

In this review paper different modeling approach and their 
comparison along theoretical and experimental results was 
reviewed. It also investigates the engine parameter and their 
effect on the performance of the engine major pollutant such as 
NO, and soot formation. Classification of various diesel engine 
models has been presented in Fig. 1. 


2. Engine submodels 


Submodels included in the study of the smallest region inside 
the engine cylinder which helps in proper understanding of the 
actual phenomena of the engine. It includes the different engine 
submodels which also help to study complex behavior of the 
combustion. 


2.1. Intake and exhaust manifold submodel 


Intake and exhaust manifold submodel deals with the mass 
flow rate inside the cylinder either the air entering during the 
intake stroke, mass of fuel during the compression stroke and the 
amount of gases flow through the exhaust stroke. The timing and 
rate of fuel injection into the chamber affects the spray dynamics 
and combustion characteristics. If the upstream pressure of the 
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Fig. 1. Classification of diesel engine model. 


Table 1 
Mass flow rate equation for flow through exhaust valve. 


SI. no. Mass flow rate equation Critical condition 


1. d 1/k-1) P, kk) ne 
as Oy HeFe Fa (= TT 24) PS KET: Supercritical flow 


k/(k—1 eas 
2. dme _ 1 p [ok (PK _pikernmy Pe > am.. Subcritical flow 
dg Sn HeFe RT \/ K-T | pak — pero ke” 


injector nozzle is known and assumed that flow through the 
nozzle is incompressible, one dimensional and quasi-steady, then 
the mass flow rate of fuel injected inside the cylinder through the 
nozzle is given by [2] 


m = CgAn,/2p, AP (1) 

In this type of mass flow rate, the intake and exhaust tanks 
(submodel) are treated as plenums with known temperatures and 
pressures. Earlier these models do not provide any information 
whether the flow is subcritical or supercritical at the inlet or 
outlet. They can be solved with the help of conservation of mass 
equation: in order to calculate the mass flow rate the equation is 
given by [15] 


g 
m=AeP RT, (2) 


2y [P PC +D 
=y | (Zar) (in 6) 


These equations are further modified by taking into account 
the critical behavior of the fluid at intake manifold and exhaust 
manifold. The intake process is assumed as a subcritical flow and 
the intake flow rate is given as [16] 


ptks +1)/ks 

rome (4) 
During the exhaust process, the flow at outlet valve is considered 

as a supercritical flow due to higher pressure difference between 

cylinder and exhaust port. When the pressure at the outlet valve 


decreases, then the behavior of flow changes to subcritical flow. 
The calculation equation is given in Table 1. 


dm 1 2ks (P/s 
= gz hFs 
dọ 6n ks—1 \ p, 


2.2. Combustion submodel 


Combustion model determines the rate of combustion. White- 
house and Way [17] proposed a semi-empirical combustion 
model for the calculation of rate of combustion. The fuel is 
injected inside the cylinder in the liquid form with conical jets 
and forming the burning zone. However, before the fuel has to be 
burned it should be first heated, evaporated as well as mixed with 
sufficient amount of air entrained from the air zone through 
diffusion. This prepared fuel is ready for the chemical process and 
may be burned according to the rate governed by a chemical 
kinetics equation [18]. For the calculation of preparation rate the 
equation used is given as 


P, = Km} Mh yp Pos (5) 


This expression was based on the assumption that the rate of 
preparation is proportional to the surface area of all droplets 
having uniform diameter, while Spalding [19] suggested that the 
rate of reaction is proportional to the droplets non-uniform 
diameter rather than uniform diameter. In order to compensate 
the non-uniform droplets diameter distribution exponent makes 
allowances for that fact. The last term in Eq. (5) shows the effect 
of oxygen availability on mixing to find out the reaction rate, for 
that purpose Arrhenius expression is used. 


KPo 4 
R= etat | P—R)d 6 
N/T i (P—R)do (6) 

The combustion is further divided into two parts as diffusion 
and premixed combustion. 


2.2.1. Diffusion and premixed combustion 

Diesel combustion is mainly considered as diffusion combustion 
however it consist of premixed and diffusion phase, controlled by 
a fuel-air mixing process. Premixed phase of combustion occurs 
before the start of diffusion combustion. It plays a significant role 
in combustion process. Premixed combustion greatly influenced 
by the engine speed, load and injection process, mainly the 
quantity of fuel injected in cylinder during the ignition delay. 
The diffusion phase is governed by the turbulent mixing of air and 
fuel. The in-cylinder turbulence resulting from air motion and 
also fuel jet induced turbulence, which is predominant in modern 
high pressure injection system is used in modern engine [20-22]. 

In the recent years different burnt rate approaches are devel- 
oped for premixed as well as diffusion combustion. Firstly, Vibe 
formula [23] was proposed for gasoline engines with totally 
premixed cylinder charge. This Vibe formula is also applicable 
for diesel engine and widely used for the diesel engine combus- 
tion modeling. This formula has only few drawbacks that it do not 
taking into account the consideration of injection rate. Another 
empirical approach was proposed by Schreiner [24] that includes 
the entire heat release history for combustion from polygonal 
sections followed by a hyperbolic function for the burn-out phase. 

Austen and Lyn [25] described the individual burn rate as 
mixing controlled process by dividing the injected fuel in a 
number of packages. This work is further extended by Constien 
[26] who used a multizone approach for droplet evaporation, 
formation and mixing and explained it in detail. The burn rate 
formulation proposed by Hohlbaum [27] using the polygonal- 
hyperbolic burn rate history was further extended by a spray 
submodel based on the fact of gas jet theory. Weisser and 
Boulouchos [28] firstly described the diffusion combustion as a 
single stage mixing controlled process. Barba [29] carried out this 
formulation by adding a premixed combustion part and also 
considered the pilot injection. 
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2.3. Fuel spray submodel 


2.3.1. Break up time calculation 

For the calculation of break-up time, Hiroyasu equation [13] 
has been most commonly used. In the quasi-dimensional com- 
bustion model of diesel engine modeling which is based on the 
theory of phase-divided spray mixing model, it is considered that 
the fuel droplets subjected to atomization period, gasification 
period and eventually become gaseous after completion of break 
up time. The break up time is calculated according to Im and Huh 
[15] is given as 


pido 
V Pa AP 


Reitz and Diwakar [30,31] included droplet break up for dense 
spray model calculations. They used KIVA-2 software code for two 
breakup regimes: (i) bag break up (ii) stripping break up. They 
shows that the droplet break up has strong influence of the spray 
penetration, vaporization and mixing in high pressure sprays. The 
size of spray droplets are controlled by coalescence and break up. 
They also include that droplet break up dominates in hollow-cone 
sprays because coalescences are minimized by expanding the 
spray geometry. Reitz and Diwakar [31] concluded that injected 
droplets having a diameter equal to the nozzle exit diameter. 
For the validation of computed length, droplet size, spray penetration 
and droplet velocity data compared with the available measured data 
in high pressure spray. 


tp = 28.65 


(7) 


2.3.2. Droplet diameter after break up 

In theory of combustion it is assumed that droplet after break 
up having the same initial diameter which is equal to the Sauter 
Mean Diameter (SMD) having neglecting droplets size distribu- 
tion and the details of the atomization process. The following 
Eqs. (8)-(10) which are commonly used for the calculation of 
SMD in multizone models are given by Hiroyasu et al. [32] to 
compensate the variation of droplet size from one fuel parcel to 
the next allowance is made depending on the operating condition 
[33]. 


as 0.54 0.18 

o = 4.12Re?-!2we,-975 (H) (Ein) (8) 
a a 

dis 0.25.) — up?) (PF 

a = 0.38Re?*> we; 932 11937 } \ p=047 (9) 

dzz d;>) dp 

dr = max (7 dh T 


Ingebo [34] developed a new theory for liquid jet breakups and 
shows that water jets in swirling and non-swirling flows, based 
on the acceleration waves. A correlation was developed for the 
ratio between the orifice diameter and SMD, and it can be 
expressed as 
dinj 
SMD 

Eq. (11) obtained is compared with the results of aerody- 
namics theory of liquid jet break up proposed by Reitz and Bracco 
[35] which predicts that SMD of liquid droplet is given by 


= C(W_Re)'/* a1) 


SMD=Bx 2 a2) 
pus; x 


inj 


2.3.3. Fuel spray dynamics 
In case of multidimensional models entire combustion cham- 
ber is taken as the computational domain and the fuel spray 


is divided into a number of zones only for multizone model. 
Each zone control volume is treated as an open system where 
mass and energy equations are solved for individual zone. 

In order to avoid the computational efficiency in case of 
multidimensional model, multizone model describe better the 
spray penetration with the help of empirical correlations instead 
of solving the full momentum equation. By the implementation of 
this method quasi-dimensional model provides the fastest and 
non-expensive means of generating the spatial information required 
to predict pollutant emission. 

Since the multizone model based on empirical correlations to 
study the spray evolution, the fidelity of the spray penetration 
model is extremely important for accuracy. For accurate predic- 
tion of spray penetration Hiroyasu and Arai [36] proposed the 
following correlations before and after break up. Hiroyasu and 
Arai [36] correlation used is listed as follows: 


S=0.39 2Ap O<t<tp (13) 
€ 
Ap 

S=2.95 p, ¥ dot th <t (14) 
a 


These correlations have been widely implemented in multi- 
zone models by different researchers [37-40]. These equations 
are valid only if the fuel elements are on the spray centerline. 
If the fuel elements in the radial direction away from the centerline 
then the equation changes to 


Xr = X exp [—4.403 x 1073(L—1)”] (15) 


where L (=2, 3, 4, 5,...), denotes the number of fuel elements 
injected per injection interval. 


2.4. Ignition delay period 


Ignition delay period is defined as the time interval between 
the actual dynamic injection point and ignition. The ignition delay 
period correlated to the cylinder pressure and temperature 
through the following semi-empirical equation given as [41]. 


ta = a, pe! (16) 


Ikegami et al. [42] and Tsao et al. [43] determined the 
induction period of diesel sprays in a rapid compression machine 
and reciprocating engine. Gerrish and Ayer [44], Spadaccini [45], 
Walsh and Cheng [46] and Scharnweber and Hoppie [47] con- 
cluded that ignition delay period decreases by the preheating of 
fuel before injection. Parker et al. [48] studied the effects of 
preheating of both the fuel and compressed air in a constant 
volume bomb to supercritical temperatures of pure n-dodecane 
which shows that ignition delay period can be expressed as 


t4 =A x e('a/T) (17) 


Parker et al. [48] shows that preheating of the fuel before 
injection has a significant role in shortening of ignition delay 
period by diminishing the vaporizing and mixing times and 
concluded that induction period is insensitive to the combustion 
chamber pressure. 

Hardenberg and Hase [49] proposed an empirical formula for 
the direct injection diesel engines to predict accurately ignition 
delay as a function of the fuel characteristics like cetane number, 
cylinder pressure, and temperature. Dent and Mehta [50] shows 
that result of this formula carried out good agreement with 
experimental data over a wide range of engine conditions. 
However, the pressure and temperature used in their correlation 
are identified as the corresponding conditions, when the ignition 
delay period starts at the end of TDC and can be estimated by 
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Table 2 
Heat transfer coefficient correlation given by different authors [2]. 


Table 3 
Calculation of rate constant [74]. 


S. No. Author name Heat transfer coefficient equations 
Eichelberg Nee =7.8 x 107!" x up333 x p95 x 795 
2. Hohenberg hen = C1 x VS x p°’ x T24 x (up?333 +c2)®® 
cy =13 x 107, & =1.4 
3. Woschni New = 3.26 x 1073 x D~? x pS x T7255 x yp8 


using a polytropic model for the compression stroke. Belardini 
et al. [51] included that ignition delay time is a function of the 
compression ratio, rpm and swirl in direct injection diesel 
engines. It has also been concluded that there is no such correla- 
tion that fits experimental data. More fundamental equation 
proposed by Hardenberg and Hase [49] to calculate the ignition 
delay period is written as 


t= 


(03640235) exp| Ea( =" + (22 A a8 
NSE See Ler | tal RT 17190) * (P-124 


Considering the effect of change of pressure and temperature 
in combustion zone, the following condition is specified [52]: 


til 
t= f ido, tę>1 (19) 


Wu [53] studied the effect of engine speed on the ignition 
delay period and concluded that as the engine speed is increased 
fuel burning rate in the premixed combustion phase decreased 
which cause increases in the ignition delay period. Wu also added 
that there is no correlation available to find the direct relation 
between ignition delay and premixed burning rate. 


2.5. Heat transfer submodels 


The heat transfer correlations of Annand [54] (including both 
the convective and radiations terms) Hohenberg, Woschni, Nussle 
and Eichelberg have been most commonly used for the calcula- 
tion of instantaneous average heat transfer coefficients. Keribar 
and Morel [55] proposed exceptional model for the heat transfer 
calculation which includes the heat transfer through cylinder wall 
as well taking into consideration the heat transfer through 
cylinder head, piston crown, cylinder liner and valves. From the 
transient heat flux measurement in diesel engine it is concluded 
that due to spatial non-uniformities of the fluid flow and 
combustion cause different temperature histories at different 
point in the cylinder head [56]. 

Further study has been carried out for the heat transfer 
through engine to coolant and found that 20-35% of the fuel 
energy is transfered to coolant [57]. Primary heat transfer 
mechanism takes place through convection from the cylinder 
gases to the surrounding areas. Also heat transfer due to radiation 
is taken into consideration because temperature inside the 
cylinder during combustion is very much higher. Heat transfer 
through convection is given by Newton law of cooling: 


e = heAc(Teyi—T wan) (20) 
The convection heat transfer coefficient (h,) can be calculated 
by using different empirical approaches such as Woschni, 
Annand, Hohenberg and Eichelberg. Annand’s formula has been 
found to be more fundamental equation which can be used to 
calculate the rate of heat transfer. This equation seems to be more 
fundamental than other alternative formulae available in the 


Reaction K ir (cm?/mol s) Author 

1 1.8 x 104 x exp(-38400/T) Baulch et al. (1991) 
0.544 x 104 x T°! x exp(-38.020/T) GRI-MECH 3.0 (2000) 
0.76 x 10* x exp(-38-000/T) Heywood [2] 

2. 6.4 x 10° x exp(-3150/7) Baulch et al.(1991) 
9.0 x 10° x exp(-2280/7) GRI-MECH 3.0 (2000) 
1.48 x 108 x T! x exp(-?860/T) Pattas (1973) 

3. 6.4 x 10% Baulch et al.(1991) 
3.36 x 10” x exp(-195/1) GRI-MECH 3.0 (2000) 
4.1 x 10” Heywood [2] 


literature. This equation considers net heat transfer as the summation 
of both radiative and convective heat transfer rate. The net heat 
transfer is given by [58] 


dQ dg ia 
a= a (Re ) (Tw-Tz) +¢(Tw-Tz) (21) 

The total heat transfer in the two zones at each angle step is 
distributed in proportion to their kilo moles and absolute tem- 
peratures [59]. The heat transfer coefficient correlations given by 
different authors have been given in Table 2. 


2.6. Pollutant formation 


Basically nitrogen oxide and soot formation are two major 
pollutants from diesel engine emissions. Nitrogen oxide and 
particulate matter emissions from diesel engines are increasing 
rapidly due to usage of high speed engines. Due to the advantage 
of direct injection diesel engine leads research to minimize the 
emission aspect of high speed engines. The main challenging job 
for the designing of diesel engine is the reduction of NOx and soot 
simultaneously. So, a lot of research has been carried out on 
design of combustion and fuel injection system to minimize the 
pollutant emission. In order to control these emissions proper 
combination is required between combustion chamber config- 
urations, fuel injection system nozzle diameter and optimum fuel 
injection timing. 

For the objective of calculating the mass of combustion 
products of each zone, the complete chemical equilibrium system 
proposed by Vickland et al. [60] is used. It is considered that at 
any instant of time for each zone 11 species are present. These 
species are O2, N2, CO2, H20, H, H2, N, NO, O, OH, CO. By using this 
scheme a set of 11 equations is formed, by which the concentra- 
tion of each zone’s (11 species) has been obtained. By solving 
these equations, mass of fuel burned, volume temperature and 
mass of air entrained in each zone has been obtained [61]. 
Chemical equilibrium based upon the assumption, so it cannot 
correctly predict the concentration of NO and this drawback is 
overcome by chemical kinetics formation scheme proposed by 
Lavoie et al. [62]. 

Besides that a lot of work has been carried out on the property 
of diesel fuel to see the effect on combustion and exhaust 
emission of variable fuel property. To reduce the effect of acid 
rain, sulfur content in the diesel fuel has been reduced. Other fuel 
properties which are related to the improvement of engine 
performance and emissions include cetane number, density of 
fuel, viscosity, distillation temperature, oxygen and aromatic 
content has also been studied by different researchers [63-68]. 

Nakatani et al. [69] used the diesel particulate NO, reduction 
(DPNR) catalyst and found that reduction in the NOx and 
particulates matter simultaneously. It has been noticed that more 
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than 80% of NOx and PM has been reduced. Hosoya et al. [70] 
proposed model with selective catalytic convertor and diesel 
particulate filter (DPF) to reduce the NO, and particulate matter 
and reported that by using selective catalytic reduction (SCR) NOx 
emissions reduced by 72%. The above used model with the 
combination of DPF has been used and found a significant change 
in the reduction of particulate matter upto 93%. 

NOx is strongly affected by higher temperature inside the 
combustion chamber during the premixed and diffusion combus- 
tion. By lowering the combustion temperature inside the cylinder 
leads to reduction of NO, emission but in other hand lowering 
temperature causes reduction of thermal efficiency which is not 
desirable. Supercharging strongly affect the combustion and 
reduces the injection rate which results in reduction of NO, 
emission. Uchida et al. [71] studied the effect of supercharging 
and concluded that favorable result are obtained in the terms of 
fuel consumption and exhaust emission. 

The extension of Zeldovich’s model with the inclusion of OH 
[72] scheme is generally accepted as a very successful method for 
predicting NO concentration. The value of temperature, volume 
and equilibrium concentration for each zone is known and by 
using these values the concentration of NO in each zone has been 
calculated by integrating the equation referring to the change of 
NO concentration at any instant of time. By the addition of NO 
concentration for all zones at a certain time, the total [NO] 
concentration inside the combustion chamber [73] has been 
obtained. Chemical reactions are given below: 


kı 


0+N2@NO+N (22) 
N-+0)¢3NO+N (23) 
N+OH<3NO+H (24) 


Rate of reaction kı, k2, k3 are calculated with the help of 
Table 3. 


2.7. Soot formation submodel 


Soot is basically produced when nondiluted spray produces 
fuel rich vapor regions and undergoes chemical reaction. It has 
also been formed on liquid film along the cylinder walls. Flower 
[75] shows that with the increase in pressure volumetric rate of 
soot formation in laminar diffusion flames increases. Rah et al. 
[76] proposed a correlation of NO and soot formation with the 
ignition of liquid fuel droplet. It also included that soot formation 
can be suppressed by delaying the ignition, vaporizing as much as 
fuel possible prior to the ignition. Rah et al. [76] also proposed a 
possible method of controlling the NO, and soot emission by 
using of low oxygen concentration in oxidation to delay ignition 
and high fuel-air ratio to reduce NO,. 

Uyehara [77] shows that in diesel engine maximum amount of 
the soot occurs in the temperature range of 2000-4000 K. It also 
shows that less than this temperature, soot formation rate is 
negligible. Kagami et al. [78] studies the effect of fuel Cetane 
number, viscosity and distillation range on the performance of 
gaseous and smoke emission, knocking, fuel consumption, start- 
ability and cold performance of diesel engine. The result shows 
that the fuel with low Cetane number produces higher emission 
and also shows that the fuel consumption is high as compared to 
fuel which is having very high Cetane number. 

Tateishi et al. [79] proposed a new combustion model for 
indirect diesel injection (IDI) diesel engines and show that this 
model reduces the fuel consumption, emission and noise. Engine 
noise is due to rise in very high peak pressure at the early stage 
of combustion, when there rapid heat release rate has been 


occurred. Sinnamon and Cole [80] used sampling to study the 
pollutant emission from a prechamber stratified charge engine. 
They concluded that prechamber induces turbulence which helps 
to increases the fuel-air mixing rates and reduces the combustion 
interval. They also show that with the increment in the degree of 
stratification beyond a certain limit resulting in higher formation 
of NO and CO emission even if lean air-fuel ratio has been used. 

Anna [81] proposed a thermodynamic model to predict the 
soot formation and oxidation rate of diesel engine. This model has 
been based upon the one step chemical reaction and D? law and 
also included droplet interference and forced convection effect. It 
has been assumed that during the vaporization and combustion 
processes liquid fuel droplets maintain their spherical symmetry 
throughout. 

Hiroyasu et al. [13] proposed a model for the calculations of 
net soot formation. According to this model, there is a soot 
formation oxidation rate, which consider the available fuel mass 
and oxygen of each zone (by including its partial pressure) 
respectively. So, net soot formation rate is calculated by subtract- 
ing the soot oxidation rate from the formation one. By integrating 
each zone separately over time, the total value of soot formation 
inside the combustion chamber is calculated by adding the 
corresponding values for all the zones. Hiroyasu et al. [13] 
proposed soot formation model is implemented in multizone 
models to calculate the soot formation. The soot formation rate 
is calculated by assuming that vaporize fuel having a first-order 
reaction which is given by the equation: 


dm 

= = Myj- M; (25) 
My = Ar Mfg pS e~(Es/RT) (26) 
m;, = Ac Ms apa e- (Esc/RT) (27) 


Lipkea and DeJoode [82] modified the previous model and 
formed a new set of equation to calculate the net soot formation. 
According to this model soot formation and soot oxidation rate is 
calculated as follows: 


dmg 


T =Ay dm? pos e-(Es/RmoiT) 28) 
dm 
dE = Ae den (pe) pre Clt 29 


Therefore, the net soot formation is expressed as follows: 
dms _ dmg  dimgc 
dt dt dt 

The net soot formation rate is proposed by Khan-Hiroyasu- 


Belardini expression and soot oxidation rate given by Nagle- 
Strickland Constable expression is given as [73]: 


(30) 


dm, dmy = dmso 


dt dt dt GD 
d 

s =Ap eT) mp (32) 
dms, _ 6M- 

de = pD; MsRso (33) 


Single zone, two zone and multizone models are proposed by 
the different author as listed in Table 4. 


3. Results and discussions 


In this literature review different modeling approach is used to 
understand the behavior of diesel engine modeling and on the 


Table 4 


Comparative study of different models of CI engine. 


Author name/country and year Model used Engine spe. (B/S/CR) Parameter studied HTC/exp Submodel used Results 

rpm val. 

Ziarati/UK (1990) [83] Single zone model (215/241/11/) na Delay Period, droplet Annand/ Combustion, heat transfer Model is capable to run under heavy fuel and diesel fuel 
penetration yes and spray model Shows good agreement of cylinder pressure and 

temperature with experimental result. 

Ishida/Japan (1994) [84] Two zone model (100/105/ na ) 3500 NO, emission, pressure, na/yes Heat release model, Ignition NO, emission is reduced by retarding the time by decreasing 
temperature and ignition delay model, NOx emission the maximum combustion temperature and pressure in 
delay model expansion stroke. 

NO, reduces by using small nozzle diameter is caused by 
decreasing the combustion pressure. 

Rakopolous/Greece (1994) [85] Multi zone model (86/80/18.5/) Insulation effect on K-EModel/ Fuel spray, evaporation, ignition No remarkable change in engine efficiency has been seen 

1200,1500,3100 performance and no delay, heat transfer, chemical but increases in the exhaust enthalpy helps to improve 
exhaust emission on combustion model overall engine efficiency by using power turbine 
diesel engine Increment in the NO and soot formation occurs by the 
insulation effect that is not desirable. 

Belardini/Italy (1996) [86] 3-D model (199/95/17/)na NOx and soot formation, na/yes NOx, soot formation, ignition The model over predicts both acetylene and soot volume 
heat delay, fraction. However, it is noted that better results can be 

obtained with a careful optimization of the model constants. 

Yong/Pohang (2000) [15] Multi-zone model (123/155/17.1/) IMEP, pressure traces, Woschni/ Intake and exhaust manifold, In cylinder flow model is combined with Hiroyasu model 


Jung/Michigan (2001) [87] 


Papagiannakis/Greece (2004) [88] 


Rakopolous/Greece (2004)[58] 


Zheng/China (2005)[89] 


Quasi - 
dimensional 
multi-zone model 


Single zone model 


Two zone model 


Multi-dimensional 
model 


1200,1400,1600,1800 


(130/160/15)2100 


(85.73/82.55/17.6) 
1000-3000 


(82/88.90/na)4500 


(95/115/20) 
1200,1600 


NOx and soot emissions, yes 
fuel injection 
timing,swirl ratio 


Pressure, temperature, Assian- 
ignition delay, Heywood/ 
combustion, NO and soot yes 
formation 


Pressure, heat release, na/no 
brake specific fuel 

consumption, NOx, CO, 

HC and soot emission 

The growth of the fuel Annand/ 


spray zone, jet mixing, yes 
C-H-O system 


Pressure, temperature, 
ignition timing, natural 


na/yes 


NOx, Ignition delay, Air 
entrainment, In cylinder flow 
effect and Fuel evaporation 
model 


Ignition delay model, heat 
transfer model, air entrainment 
model, spray evaporation model 


HCCI model, heat transfer 
model, pollutant formation 
model 


Fuel Preparation,combustion 
Model, NOx, Soot Formation 
model 


both model shows results consistent with the general 
operating characteristics of a DI diesel engine. 

The calculated imep, pressure traces, brake specific NOx and 
soot emissions studied by the present model shows good 
agreement with the test bed measurements at full and part 
load conditions for a range of the engine speeds. 


Extended Zeldovich mechanism accurately predicted the NO 
and soot formation. At least five zones are recommended to 
determine the NO and soot in radial direction. 

Modified Hiroyasu correlation for spray tip penetration 
show good agreement with experimental data in a range of 
injection pressure. 


Exhaust emission and performance of engine study have 
been carried out at different speeds and loads for diesel and 
dual fuel (combination of diesel and gaseous fuel). 

Dual fuel combustion process has lower peak cylinder 
temperature in comparison to conventional diesel. 

NOx and soot emission is low and HC and CO emission is 
high in dual fuel operation compared to normal diesel 
operation. 


Two zone models proves to very efficient for the 
determination of engine performance and exhaust emission 
by concerning the effect of two major parameter of load and 
injection timing. 

Theoretical result show good agreement with exp. value. 


Injection time, natural gas composition and initial 
temperature have a significant effect on ignition, 
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Table 4 (continued ) 


Author name/country and year Model used Engine spe. (B/S/CR) Parameter studied HTC/exp Submodel used Results 
rpm val. 
gas composition, NO Combustion model, turbulence combustion and emission. 
emission model, mixture formation Increased amount of ethane in natural gas provide advance 
model, No formation model in ignition and increment in NO emission. 
Auto-ignition of natural gas occurs at 1300 K of mixture 
temperature and NO emission reduce when fuel injection 
timing delay. 
Tauzia/France (2006)[90] Phenomenological (123.8/165/16.5) Premixed combustion, na/yes/cfd- Premixed combustion model, Model studied at injection timing -16° to +12°ATDC shows 
Model 1900-3000 pressure, temperature kiva 3v Ignition delay model good agreement with experimental and CFD analysis from a 
quantitative and qualitative point of view for ignition delay 
and in-cylinder pressure evolution during combustion. 
Xingcai/Chima (2006) [91] Multi-zone model (98/105/18.5)1800 IMEP, thermal efficiency, na/yes HCCI combustion model, Energy input, ignition timing, maximum gas pressure is 
pressure, heat release pollutant model deteriorated with increase of ethanol addition. 
rate, ignition timing, For all fuels, CO emission is very high at 1.5-2.5 bar of 
emission characteristics engine load. 
When IMEP is larger than 3 bar ignition timing lies in -5°CA 
ATDC-TDC. 
Chmela/Austria (2007) [92] Zero Dimensional na Premixed combustion na/no Ignition delay,Premixed and An improved model for diffusion combustion is presented 
ROHR Model based on Arrhenius and diffusion combustion,Reaction which shows a sequential approach. 
Mangussen equation. Rate combustion Model Ignition delay and premixed combustion model solved by 
Arrhenius and Magnussen approach shows better 
description of the effects of considerably transient ROI 
shapes by using a more detailed consideration of the 
processes in the fuel jet. 
Sundarapandian/India (2007) [93] Multizone model = (87.5/110/17)1500 Heat release, work done, Annand/ Combustion and heat release, Heat release is reduced by 4-8% for vegetable oil esters 


Machrafi/France (2008) [94] 


Sahin/Turkey (2008) [95] 


Single zone model 


Multi-zone model 


(82.55/114.5/4-14)na 


na 


and harmful pollutant yes 
such as HC, CO, NOx and 
smoke. 


Heat release ratio, 
equivalence ratio, 
ignition delay, 

compression ratio 


na/yes 


Details of fuel spray, heat Annand/ 
transfer, swirl yes 
model,engine 

performance 


heat transfer, engine 
performance like specific fuel 
consumption, brake thermal 
efficiency 


Combustion model ignition 
delay model, heat released 
model 


Heat transfer, combustion 
model, spray model, NO 
emission model, 


compared to diesel fuel. Pressure of vegetable oils esters are 
reduced about 4%, 5% and 7% respectively for Jatropha, 
mahua and neem oil ester compared to diesel. 

But the negative impact has been seen on thermal efficiency 
reduced by 3%, 4% and 5% for Jatropha, mahua and neem oil 
ester compared to diesel. 

CO reduced by 19% for Jatropha and 16% for the mahua and 
neem oil ester compared to diesel. 


Heat release increase with increase in equivalence ratio 
which is affected by ignition delay. 

Compression ratio provides advancing in ignition delay and 
lower residual gas fraction. 

Variation of inlet temperature is directly affected the 
ignition delay but heat release distribution did not affected 
the ignition delays. 


Indicated thermal efficiency at 2500 rpm by present model 
is 0.4420 and Li model is 0.43 shows that in result only 3% of 
error occurs. 

Equivalence ratio changes between 1.086 to 1.1.49 while 
actual ration like 0.650-0.80 so model is not adequate at low 
equivalence ratio. 
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Semin/Malaysia (2008) [96] 


Garcia/Spain (2009) [97] 


Huang/China (2009)[98] 


Kannan/India (2009) [99] 


Komninos/Greece (2009)[100] 


Kuleshov/Moscow (2009) [101] 


Liu/China (2009) [102] 


Single zone model 


Single zone model 


Multi-dimensional 
CFD model, single 
zone model 


Single zone mode 


Multizone model 


Multizone Model 
(RK-model) 


Two Zone Model 


(86/70/na)200-4000 


(95/100/19) 
1500-2400 


(115/115/17)1400 


(87.5/110/17.5)1500 


(120.65/140/17)na 


(88/85/17.6)1800 


(102/120/18:1) 
800-2000 


MEP, mean piston speed, na/no 
specific power output, 

sfc intake valve mach 

index. 


Pressure, fuel specific 
consumption, start of 
combustion, NOx, CO, HC 
and soot emission 


na/no 


Pressure profiles, 
combustion process, 
HC and CO emission 


na/yes 


Heat release rate, Hohenberg/ 
pressure, brake thermal yes 
efficiency, bsfc NO 

formation 


Mass transfer effect on Annand/ 
emission i.e. unburned yes 
HC and CO emission, 


Heat release inside the | Woschni/ 
cylinder and NOx yes 
formation, inlet manifold 

pressure and EGR rate 


Mass distribution, heat |= Woschni/ 
transfer,mass exchange na 
between zones aheat 

transfer in porous 

medium and studied the 

effect of PM on 

performance 


Parametric model -MEP, SFC, 
power output, mean piston 
speed. 


Combustion model, pollutant 
model, HCCI combustion model 


Turbulence model, combustion 
model, mixture formation 
model, HCCI combustion 


Ignition delay model, heat 
transfer model, pollutant 
emission model. 


Heat transfer as well as 
pollutant formation model 


combustion model, fuel spray 
modelgas exchange model, 
ignition delay model 


Heat transfer model, mass 
exchange model, porous 
medium model, chemical 
kinematics 


From engine modeling it is concluded that highest brake 
power is 4.314 kW and indicator power is 5.725 kW occurs 
at 3000 rpm. 

The engine modeling brake torque maximum is 14.3662 

N m and indicated torque is 18.743 Nm occurs at 2800 rpm. 
At 2800 rpm, the engine model has shown the lowest value 
of the BSFC and the indicated specific fuel consumption. 


HCCI combustion process provides low NOx and soot 
emission than the diesel mode. HCCI with exhaust gas 
recirculation (EGR) has a great impact on combustion and 
emission. 

When angular speed increases, start of combustion that 
affected by inlet temperature at constant EGR rate, decrease. 


HCCI combustion process does not take place 
simultaneously in combustion chamber core zone. High 
temperature reaction occurs in vicinity of piston surface. 
Unburned fuel and CO increase when Di-methyl ether 
(DME) equivalence ratio decreases. 

Pressure profile of multi-dimensional CFD model is more 
accurate than the single zone model because turbulence 
flow and wall heat transfer is measured. 


From simulation result it is concluded that 18% and 21.5% of 
reduction in NO is achieved with 10% and 20% dilution of 
diesel with water emulsion respectively. 

The brake thermal efficiency increases as reported by Abu - 
Zaid by 3.5% by increasing the water emulsion by 20% over 
the use of diesel for the engine speed range studied. 


Neglecting mass transfer prior to main heat release would 

result in a severe underestimation of the HC accumulated in 
the crevices during combustion (60% for the case Studied). 

Neglecting mass transfer during combustion and expansion 
would result to and underestimation of the final HC amount 
by at least 14% for the case studied. 


The RK-model accounts for drop size interaction of free 
sprays with swirl effect, spray and wall impingement, 
evolution of near-wall flow formed by spray, hit of fuel on 
cylinder head surface and swirl intensity on heat 

release rate. 

This model helps in designing of piston bowl shape, 
numbers of injector holes, diameters and directions in 
match with swirl intensity to control emission and fuel 
consumption. 


The compression ratio of 13 is great enough for the auto- 
ignition in the porous medium engine, which is not so 
crucial like in the HCCI engine. 

For the same intake pressure of 1 atm, the pressure in the 
PM engine is relatively lower than that of a traditional 
engine due to the voids in PM, which lowers the ignition 
delay and higher engine output. 

Due to higher average temperature in PM engine NO 
emission is higher at excess air ratio of 1.6. 
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Author name/country and year Model used Engine spe. (B/S/CR) Parameter studied HTC/exp Submodel used Results 

rpm val. 

Dhuchakallaya/UK (2010) [103] Spray model na Combustion, spray na/na Combustion model, Ignition Higher lift of length cause the increases of air entrainment in 
dynamics and ignition model the premixed combustion zone resulting in high reaction rate. 
delay Simulation results show good agreement with experimental 

results at time of 4 ms of start of injection. 

Gogoi/India (2010) [104] Single zone model (87.5/110/14- Engine performance- Annand/no Heat transfer, ignition delay Maximum brake power occurs at design speed and further 

20)600-2000 IMPEP, BTE comparative model, gas exchange model increases of speed decrease brake power. For CR of 17.5 the 
study of diesel with maximum peak power occurred at 1500 rpm for B20 and B40 
biodiesel blends blends while for B60 it occurs at 1400 rpm. 

Increased brake power was observed for B40 and B60 blends 
as compared to diesel fuel. 

Blends B20, B40 and B60 shows an increase in thermal 
efficiency due to increases in break power as compared to 
diesel engine. 

Prasad/India (2010) [105] Two Zone Model = (11/127/16)1000- Pressure, heat release, Annand/ Preparation rate and reaction The cylinder pressure for low heat rejection biodiesel is 

1500 heat transfer and yes rate model for the calculation of lowers than LHR diesel about 1% and higher about 4% and 
performance heat transfer rate. 2.5% for TC biodiesel and TC diesel. 
characteristics -SFC, BTE Brake thermal efficiency of LHR biodiesel engine is lower 

than that of LHR diesel and higher than TC biodiesel and TC 
diesel operation this is because of lower calorific value of 
biodiesel as compared to diesel fuel. 

Qi/China (2010) [106] Quasi-dimensional (135/150/17) 1500 Diesel engine Newton/ Spray mixing, combustion Simulation result shows good accordance with experimental 


Rakopolous/Greece (2010) [107] 


Shankra/India (2010) [108] 


Shi/USA (2010) [109] 


two zone model 


Quasi-dimensional (85.73/82.55/17.61) 


Multidimensional 
model (CFD 
model) 


Multidimensional 
model 


500-2500 


(130/150/15.5:1) 
2000 


(137.16/165.1/16.1) 
1300 


performance- sfc,effect yes 
of air entrainment and 
spray mixing model 


Cylinder flow field and 
gas temperature 
distribution. 


na/CFD 
model 


CFD model proposed to 
study the effect of intake 
pressure, fuel injection 
on the performance of DI 
engine 


Eckert/yes 


Pressure, injection 
timing, specific fuel 
consumption, soot 
formation, NOx, CO, UHC 
emission 


na/yes 


model, ignition delay period, 
heat transfer, pollutant 
formation, intake and exhaust 
flow model 


heat transfer, fuel injection rate, 
spray penetration, evaporation, 
combustion and pollutant 
formation 


CFD based—combustion and 
ignitions model, turbulence 
model, NO, and soot formation 
and droplet breakup model, 


Adaptive multi-grid chemistry 
model, combustion model 


result proves the validation of phase divided spray 
mixing model. 

The relative error between simulation and experimental 
result for the specific fuel consumption and NO emission 
occur to be 2.8% and 9.1%. 

By using this model computational time is 36 s to analyze 
engine working process. 


Both quasi-dimensional and CFD model are compared for 
the different piston bowl geometry changing from (d/D) 
from 64% to 44% at different rotational speed. 

It is concluded that both models predict the quite similar 
cylinder pressure and temperature, traces whole closed 
cycle of the engine. 


It is observed that by advancing injection time cause results in 
the increases pressure, temperature heat release rate, 
cumulative heat release and NO, emissions by 6.88%. If we 
reverse the trend then NO, emissions by 6.85%. 

Higher NO, (15.03% and 58.69% at 1.21 bar and 1.71 bar 
respectively) and lower soot emission (8.82% and 51.47% at 
1.21 bar and 1.71 bar) concluded that optimum injection time 
and intake pressure is 12° bTDC and 1.21 bar. 


At high load, high volatility fuel as gasoline and E10 have 
less sensitivity to fuel reactivity. So these are better for fuel 
economy. 

For cleaner combustion, gasoline fuels are more efficient 
than normal diesel fuel. 
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Som/USA (2010) 
[110] 
na 


Reiter/USA (2011) [111] 


Boretti/USA (2012)[112] 


Komninos/Greece (2012) 
[113] 


Visakhamoorthy/Canada (2012) [114] 


Single zone model (50/200/na) 
na/yes 


Multi-zone model (106/127/17) 1000 


Multi-zone model (131/158/18)1900 


Single zone, multi- na 
zone and 
multidimensional 
model 


Multi-Zone model (72/73.6/22) 1500 


Spray and combustion 
characteristic, fuel 
penetration, 
temperature contours, 
flame index 


Pressure, heat release 
rate, engine break 
specific fuel 
consumption, soot 
formation 


Cumulative heat release, 
pressure, heat release, 
pollutant emission 


Chemical kinematics 
modeling, different 
reaction mechanisms for 
biofuel combustion 
simulation 


Pressure, heat release 
rate, Emissions, 
predicted NO, levels 


na/ no 


na/no 


na 


Yes/ yes 


KH breakup model, spray and 
combustion model, 
aerodynamically induced 
breakup model, soot formation 


Ignition delay model, 
combustion and heat release 
model, pollutant formation 
model, soot formation model 


Combustion model with jet 
ignition and diesel ignition, 
stochastic reactor model (SRM), 
soot model 


HCCI combustion model, direct 
injection HCCI, fully premixed 
HCCI, emission model 


HCCI engine model, mixture 
composition model, pollutant 
model 


Normal diesel is affected by first injection amount but 
gasoline fuel is affected by second injection timing. 


-KH-ACT model have better performance of parameters 
such as injection pressure, ambient density and fuel 
temperature with vapour penetration for evaporating 
spray. -The inclusion of turbulence and cavitations 
enhances the primary breakup process, decreases in liquid 
penetration and increases in radial dispersion. -Enhancing 
spray breakup provides ignition length and flame 
stabilization near nozzle exit. 

In dual fuel engine, ignition delay is longer due to increased 
ammonia and resulting peak combustion pressure is low. 
BSFC of diesel is high below 40% of diesel fuel energy. 

NO emission is low by using less than 40% ammonia in dual 
fuel compared 100% of diesel. 


Diesel injection during compression and expansion to build 
up the pressure are best options to develop a pressure 
profile. 

Thermal efficiency of diesel injection and hydrogen jet 
ignition combustion models is 10% larger than the diesel at 
high loads. 


Present the effect of various configurations in single zone 
model such as negative valve overlap, heavy EGR, fuel 
reforming on biofuel HCCI combustion. 

Study the multi-zone models for mixture thermal and fuel 
stratification which provide better and more realistic 
estimates on combustion duration, peak combustion 
pressure and emissions formation. 

Multi-dimensional models provide higher spatial resolution 
of the combustion chamber than the multi-zone model. 


The model represents good quality HRR and pressure 
prediction curves after calibration. 

When engine operating point is close to misfire limit, model 
is unable to handle cyclic variability i.e. to predict the 
pressure. 

Predicted NO, emission increase with small increase in 
cylinder temperature. 


HTC—heat transfer correlation, B—Bore of cylinder, Spe.—Specification, Exp.—Experimental, S—Stoke length, Val.—Validation, CR—Compression ratio. 
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basis of these modeling approach, following result has been 
obtained and discussed below: 


1. Insulation effect on the cylinder wall is studied by the multi- 
zone model and it is concluded that by the increasing the 
degree of insulation causes increment in the NO, emission and 
simultaneously increases soot formation. This increase in 
emission and soot formation is not admirable because of 
negative impact on the environment. However, the engine 
efficiency is improved by providing more insulation on the 
cylinder wall. So, it is suggested that degree of insulation 
should be properly selected. 

2. Rate of heat release (ROHR) model proves very effective to 
accurately predict the different instantaneous parameter 
which controls the diffusion combustion and also it takes very 
less computational time. 

3. Single zone approach by using the water diesel combustion 
emulsion has been successfully modeled to predict the engine 
performance parameters. It has been found out that with the 
20% increases in the water emulsion brake thermal efficiency 
(BTE) of the engine increases by 3.5% which is significant 
improvement for diesel engine. Similarly the NO, emission 
reduced from 21% to 18.5% by increasing the water emulsion 
from 10% to 20%. 

4. CFD study is carried out for the detailed analysis on the 
multizone model to study the effect of intake pressure, injec- 
tion timing on the NO, and soot formation. From the analysis it 
is concluded that with the increase in the intake pressure NO, 
emission increases while soot formation rate decreases. So, 
from the three intakes pressure range 1.01, 1.21 and 1.71 bar 
the 1.21 bar gives more desirable optimum result which 
compensate both type of emissions. Similarly out of 16, 12 
and 8 bTDC injection timing the 12 CAD found to be providing 
more optimum result of pollutant emission. 

5. Ignition delay model is studied at different cetane numbers 
and variable speed by using Hardenberg equation and con- 
cluded that with the increases in the engine speed the ignition 
delay period increase while with the increase of the cetane 
number. 


4. Conclusion 


In this present article single zone, multizone, and multi- 
dimensional models of diesel engine combustion has been 
reviewed. Single zone model has been found to more effective 
tool for immediate prediction of engine parameter but case by 
case adjustment of the wiebe function parameter is required. 
Theses wiebe function parameter deals with engine geometry 
and operating condition. Multizone model account for air 
entrainment and inhomgeneities of the mixture and overcome 
the drawback of single zone model. Multi-dimensional models 
account for both engine geometry temporal as well as spatial 
variation of the flow field. These accurately predict the spray 
behavior and fuel spray dynamics model. These model 
included the CFD analysis which accurately predict the com- 
bustion behavior and help to reduce the pollutant emission at 
great extent make diesel engine more reliable and fewer 
impact on the environment. 
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